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О задаче оптимального управления в моделях роста 
патологических образований на основе динамики Рихарда 


У даній роботі запропоновано модель керування патологічним процесом, розвиток якого описує динаміка 
Ріхарда. Модель представлено як задачу оптимального керування системою нелінійних диференціальних 
рівнянь. Сформульовано необхідні умови оптимальності та аналітично показано вигляд оптимального 
розв'язку на деякому інтервалі часу. 

Ключові слова: оптимальне керування, принцип максимуму Понтрягіна, 

необхідні умови оптимальності. 
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В данной работе предложена модель управления патологическим процессом, развитие которого описьтваєт 
динамика Рихарда. Модель представлена как задача оптимального управления системой нелинейньх 
дифференциальньх уравнений. Сформулированьт необходимьте условия оптимальности и аналитически 
показан вид оптимального решения на некотором интервале времени. 

Ключевьге слова: оптимальное управление, принцип максимума Понтрягина, 

необходимьгюе условия оптимальности. 


Вступ 


При проведенні системних медичних досліджень виникають питання прогнозу- 
вання кількісної та якісної поведінки захворювання. Це, в першу чергу, форма патоло- 
гічного процесу, на яку впливає цілий ряд невизначеностей - час формування каскаду 
специфічних плазматичних клітин, вплив ушкодженого органа на імунну відповідь, 
схема проведеного лікування та ін. Розв'язування проблем такого роду вимагає роз- 
робки відповідних алгоритмів системного аналізу |1-10). 
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Дослідження проблеми розвитку патологічних утворень у людському організмі є 
актуальним і в наш час. Зокрема, в роботах |1-7| розглядають ріст пухлинних популяцій 
на основі динаміки Гомпертца. В роботі |8| вивчають розвиток загального патологіч- 
ного утворення на основі динаміки Ріхарда. В |9| розглянуто питання стійкості в моделі 
росту патологічного утворення на основі динаміки Ріхарда. Проте не досліджене пи- 
тання керування процесом росту патологічних утворень. 

Метою роботи є розгляд задачі оптимального керування процесом росту пато- 
логічного утворення, який описує динаміка Ріхарда, та визначення необхідних умов 
оптимальності. 


Матеріали та методи дослідження 


Нехай функції п (/)  - 1, У є розв'язками системи диференціальних рівнянь 
і й 2 


щу 4з ОТ ОЗ о - о | (1) 


лч1 І 
А (г, ) - їі 
де 7, (г , а; (г) -- кусково-неперервні функції на відрізку (/,.Т), 8(г) - Мхт- 
матриця з елементами 5, (г), що є також кусково-неперервними. При цьому розв'я- 


зок п (г) розуміється в загальному, тобто: 


РЧР ДОЛОТЬЗУХ 29) а. 


З; 
ЖО, 
8 


У рівнянні (1) зробимо заміну х г )- . Тоді отримаємо: 


ах, 2 их, є Ум 05, аб (220) | о) 
0, 
Розглянемо керований об'єкт, що описується системою рівнянь (2), де х є К" - 


фазовий стан об'єкта; и є К" - керування. 


Розглянемо задачу оптимального керування 
Т 


Ли)з рю іп, (3) 


Ок и, оку ав хо яд дек, 4) 


х(0) вхо» х(Т) шх, (5) 
пе У чат, (6) 
дез и(г) -- вважаються кусково-неперервними функціями на |/,, |, 1 Точки Хо, 
х; - Задані, / є Е" - не залежить від часу і фазові обмеження при /, 2/1 2 Т відсутні. 
Необхідно знайти таке припустиме керування ий), що переводить систему з 
фазового стану х(0)- х09 у фазовий стан х(Т )- хФ , причому відповідний припусти- 
мий процес (х(г), и) надає мінімального значення функціоналу (3), де функція и'(І) - 
неперервна за сукупністю усіх змінних. 
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Застосуємо принцип максимуму Понтрягіна |10)| до задачі оптимального керу- 
вання (3) - (6). 

Введемо допоміжні змінні (г) - (цу (г),., (г) є Е" і сталу ш, - -1. Визначимо 
функцію Гамільтона-Понтрягіна: 


НСхО)и(т),г) з У Г'(хИ), нау) ни ГОСТ) (1) з 
120 (7) 
зон) у (пх, 95. Ь, (ги, (г)чк 3 а, «і -х, (г)) 


Парі (и(г), х(1)), г, 212Т поставимо у відповідність систему диференціальних 
рівнянь 


ау, | дн(х() ит) | од (хИ), пт)з1) 
ФО дх Й 2, хо 1 


І І 


із0,п. (8) 
Продиференціюємо функцію Гамільтона-Понтрягіна (7) за змінною и(1): 


5; М 


ГЕТТО; 


З умови 
з ди 


- 0 знайдемо: 


цу я зп (Ок (у 9) 


Продиференціюємо функцію Гамільтона-Понтрягіна (7) за змінними м, і х,: 


нин) нд 


аї 


- ,із0,М, (10) 
-фх 25, ь (Би уч Уа (61-х, ої 


дно С Р У 
дх Ро д|-Я бх РМ 


ізі і 


, ,із0,м. (1) 
сип но, ь (Би уч Уа (01-х, ої -пх(бууУ а (1) 


Тепер співвідношення (10) і (11) з урахуванням и(/) можна переписати у ви- 
гляді гамільтонової системи: 


т 


пу А 227 ом, (ох (бупуч Уа б -х, ої «б, 


- -- ТЕ пуд будує Уа Ух, ої - (12) 


т 


тору а, (у із уп. 


Отже, можемо сформулювати необхідні умови оптимальності для задачі (3) - (6). 
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-. Ж Ж 5 ; 
Теорема. Нехай (хХН (г),и (ї)) - розв'язок задачі оптимального керуван- 


Х ж ж - - . 
ня (3) - (6). Тоді (х (г), и (г)) є розв'язком крайової задачі (12), (5). 
Далі розглянемо скалярний стаціонарний випадок, увівши позначення 


М М 
Уа )з а У Ь ()- Ь. 
71 іч 


Припустимо, що а,р»д. 
Розглянемо систему: 


1 
е. пах(1) - пах () з --п'Ь'х (ши) 
аї 2 13) 
Чи по 1 ро 2 
і -пац(т) - 2пах(Т)и (1) - 5 п'Б'х(уг (1) 
Для визначення двох сталих із загального розв'язку системи (13) маємо дві умови (3). 
Отримали крайову задачу (13), (5) - для системи двох звичайних диференціаль- 
них рівнянь. 
Система (13) має загальні розв'язки: 


Р, о Ко БОНН (3 ПУ і А (14) 
дасте з ус! 
є о віє Дп'а С) зе НО) Др'а'еОО ра" і 
- 2пах(1) - Зпах (І) УЖНУ ці 
(з 222.2 . 
п'Б'х"(Г) 
ас? 
р а|хи)- і дай 
С, )С 
а ) РА б ЕВ оїнс, З рАп'а"!- тессі аб) 


- 2пах (І) - 2пах ") я хо х(Т)) 


се пвх) 
Розв'язок Р, - тривіальний. 
Розглянемо розв'язок Р. Підставимо крайові умови (5) в загальний розв'я- 
зок (15) і визначимо сталі інтегрування С, та С,: 
С 22,де 2 - корінь полінома 
О (7) «(4хуВ'АхІВ'-Ахх В АВ 
нОокухі В Ах Вк око у (В - ап) «0, 
а Ю-е?,де р- корінь квазіполінома: 
ОР) а -АхЦе? У рі Ах?) р нах рі нах (є?) різ 
З хіхоТ ан (е? ур» - хіх ап Че? У рі - Дух Тапе? ) ре з 


З Зхух Та? (є?) р? -4хухіТ'а'п'є рІ- Ахіх Тан (е?) ра З 
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вхо Т'А?нЧе?У р» -Ахух Т'а'п'є? р' з хух Та"'п'р'- 
4 З 
-хохТ'а" пе?) З 0 туя -бхух Та Че?) З 
ж Ахох Т'а"п'є? - хх Т"а"п 
Коефіцієнт С, визначається через значення С, і дорівнює: 


с.аці118 та? «4032 544 - 2х а С С зх Сіпіа' ЄВ б 


2 5 1: 
2 меч 


Визначивши С, та С, ми отримаємо частковий розв'язок задачі (13), (5). Вираз 
для х(Г) і Щ(Т) з урахуванням значень С, та С, підставимо у (9) і отримаємо керуван- 
ня и(1). 

Розглянемо розв'язок Р, . Підставимо крайові умови (5) в загальний розв'язок (16) 
і визначимо сталі інтегрування С, та С,: 

С, «2 ,де 2 - корінь полінома: 
О, (7) з (Ах В' Аху В" хх В Ах В) хх ох В хх ВУ 
нах Воді В-хух 2 (Ю -Тап об, 
а Вае?,де Р - корінь квазіполінома: 
ор ах ур у4х(е ри  4ххуедр' -Ахух(е р щ 
хіт Т'ачп'(е? і да дхіхТ а о бліх Т'а'пе? ух 
-АхухіТ"а"п'є? в хохіТа"п'- хохіТ а"п Чер у різ зх ГЗа'яЧе? ) ра З 
жаху Т'а"п Че? ур вух Та? (е? У р? жаху Т'а"п'є? різь 
З хіх ТапЧе?) ра -8хіх Га? (е?У з дхіх Т'а?п(е?)- р кана), 
Коефіцієнт С, визначається через значення С, і дорівнює: 


І - 4Ахупда «4СІ 4 2хупда 2 з С хіпда? С2 


Сі ат й 
и з; (Ах п'а" - Ахт'а С2зв'ї» 


С: 


Визначивши С, та С, ми отримаємо частковий розв'язок задачі (13), (5). Вираз для 


х(г) і шм/(Г) з урахуванням значень С, та С, підставимо у (9) і отримаємо керування (1). 


Висновки 


Сформульовано задачу оптимального керування. Знайдено загальні розв'язки 
крайової задачі, яка описує процес росту патологічного утворення на основі динаміки 
Ріхарда. Застосовано принцип максимуму Понтрягіна до задачі оптимального керування 
та отримано необхідну умову оптимальності. 
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арріїсайоп5 5исі а8 Біоіодіса! 5узіетл8, сопаттипісайоп пебуогк5 апа 5осіо-есопотіс зузіетя 
еїс. Аз8 а гезиії, плоге апа тоге реоріе млі! Бепебі стеайу Бу Іеагпіпе о 50Їуе Бе орійтаї 
сопіго! ргобіетя питегісаПу. 

Ттеайтепі ої а райпоєепіс Фізеа5е ргосе55 15 іпіегргеїед а5 бе оріта! сопіго ої а 
дупатіс 5узіет. Емоїшйоп ої ре дї5еа5е 15 спагасіегі7ей Бу а поп-Ппеаг огдїпагу ФіНегепнаї 
едпайоп Шаї дезсгібеє сопсепітайоп5 ої раїпогеп5, ріазта сеїя, апа апіїбоаіе5, аз ме а8 а 
питегіса! іпфісайоп ої райепі реа. М/Фоці сопіго!, Ше дупатіс плодеї еуідепсез 58иі- 
сПппіса! ог сійпіса! десау, сргопіс 5:арійгайоп, ог ппгезітаїпеа Їеїра! ого ої їре рафогеп. 

Тре дупатіс едпайопя аге сопігоПед Бу Шегареціїс асепія Шаг аНесі Ше аїе ої срапее 
ої зузіега уагіабіе5. Орбіта! сопіто! 50Пийїоп5 Краї деїеаї їБе раїроєеп апа ргезегуе ограп 
реаїйі аге депопяігаїед. Пі і5 5Ппомуп Має сопіго! ШФеогу сап роїпі Ше ухау їОууага пеуу 
ргоїосоїз їог ігеайтепі апа геплейтайоп ої ритап 4і5еа5е5. 

Іліегасйоп ої Іа5зег птайіайоп апа Біоіовіса! зузіет5 аге іпуе5бдаїед. ТПе раїроіодісаї 
Юогтабйоп огом'Ю Базед оп Вісрага5 дупатісз аге соп5ідегей. 

Ропітуавіп плахітит ргіпсіріе 15 аррпед їп Пі5 ууогК р Ппа Фе оріїта! сопіто! об 
Фе рафоіоєсісаї Гогптпайоп єгом/'їр. ТрПе єбепега! 50Гийоп ої Бе сопігоПед зубіет 15 
обіаїпей. ЕїпаПу, питегіса! ехатріе аге ргезепіед іо ехріаїп Фе апаЇуйсаї 5оїийоп. 


Стаття надійшла до редакції 01.10.2012. 
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